

  ss = (/ "ANN","DJF","JJA","MAM","SON" /) 

  sitenama_list = asciiread("radon_sitename.dat",-1,"string")   

  ns = dimsizes(sitenama_list) 


  str = new(10,"string") 
  str = "" 

  do it = 0,5-1 
  do is = 0,ns-1 

  sitename = sitenama_list(is) 
  sitenama = sitenama_list(is)  
  season = ss(it)  
  siteid = is
  passid = 0
  H1215  = False

  str(0) = "iseason = " + it
  str(1) = "isite   = " + is
  str(2) = "sitename = sname(isite)" 
  str(3) = "sitenama = sname(isite)" 
  str(4) = "season = ss(iseason)   " 
  str(5) = "siteid = isite         " 
  str(6) = "passid = 0             " 
  str(7) = "H1215  = False         " 

  asciiwrite("load_site.ncl",str) 

  system("ncl plotbox.ncl") 

  end do  
  end do 

